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Abstract 

The nucleon mean-field potential has been thoroughly investigated in an extended Hartree-Fock 
(HF) calculation of nuclear matter (NM) using the CDM3Y3 and CDM3Y6 density dependent ver¬ 
sions of the M3Y interaction. The single-particle (s/p) energies of nucleons in NM are determined 
according to the Hugenholtz-van Hove theorem, which gives rise naturally to a rearrangement 
term (RT) of the s/p potential at the Fermi momentum. Using the RT obtained exactly at the 
different NM densities and neutron-proton asymmetries, a consistent method is suggested to take 
into account effectively the momentum dependence of the RT of the s/p potential within the stan¬ 
dard HF scheme. To obtain a realistic momentum dependence of the nucleon optical potential 
(OP), the high-momentum part of the s/p potential was accurately readjusted to reproduce the 
observed energy dependence of the nucleon OP over a wide range of energies. The impact of the 
RT and momentum dependence of the s/p potential on the density dependence of the nuclear 
symmetry energy and nucleon effective mass has been studied in details. The high-momentum tail 
of the s/p potential was found to have a sizable effect on the slope of the symmetry energy and 
the neutron-proton effective mass splitting at supranuclear densities of the NM. Based on a local 
density approximation, the folding model of the nucleon OP of finite nuclei has been extended to 
take into account consistently the RT and momentum dependence of the nucleon OP in the same 
mean-field manner, and successfully applied to study the elastic neutron scattering on the lead 
target at the energies around the Fermi energy. 
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I. INTRODUCTION 


The nucleon mean-field potential or single-particle (s/p) potential is the most important 
quantity determining the s/p properties of neutrons and protons in the nuclear matter 


(NM) as well as in the finite nuclei 


0 , 


and it has been on the focus of recent many- 

the NM 


body studies of the NM, like the Brueckner-Hartree-Fock (BHF) calculations o 
starting from a realistic choice of the free nucleon-nucleon (NN) interaction 


or the 


mean-field studies of the NM on the Hartree- 


ock (HF) level, using the different versions of 


the effective (in-medium) NN interaction 5-7]. Quite interesting aspect of the s/p potential 


is its direct connection with the nuclear symmetry energy 


( 4 , 7], a key quantity necessary 


for the determination of the equation of state (EOS) of the asymmetric NM |8|. Many 
microscopic studies of the EOS were done based either on the nonrelativistic or relativistic 
mean-field potential given by realistic two-body and three-body NN forces or interaction 
Lagrangians |q, 10]. Such microscopic many-body studies did show the important role 
played by the Pauli blocking effects as well as the increasing strength of the higher-order 
NN correlations at the high NM densities. These medium effects are usually considered 
as the main physics origin of an explicit density dependence embedded in the different 
versions of the effective NN interaction, being used currently in the HF calculations of the 
nuclear structure or nuclear reaction studies. Among them, quite popular are the density 
dependent versions of the M3Y interaction (originally constructed to reproduce the G-matrix 
elements of the Reid [ll] and Paris 12] NN potentials in an oscillator basis), which have 
been successfully used in the HF studies of the NM 13Ml7| as well as in the folding model 


studies of the nucleon-nucleus and nucleus-nucleus scattering 


18 


23]. 


With the phenomenological density dependence of the M3Y interaction parametrized to 


give a realistic description of the NM saturation properties 


HQQ 


within the HF frame, 


the HF nucleon optical potential (or the high-momentum part of the s/p potential) was 


used 


13|, 


to determine the explicit energy dependence of the density dependent M3Y interaction 
24] based on the observed energy dependence of the nucleon optical potential (OP). 


These density- and energy dependent versions of the M3Y interaction have been further 


used in the folding model calculation of the nucleon-nucleus and nucleus-nucleus OP 


22], The simple assumption for the s/p potential in the NM made in Refs. 


13 


equivalent to the microscopic s/p potential of the Brueckner-Bethe theory 25], which lacks 


n 


24] is roughly 
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the so-called rearrangement term that arises naturally in the Landau theory for the infinite 


Fermi systems 


26j. Such a rearrangement term (RT) also appears when the s/p potential is 


evaluated from the total NM energy using the Hugenholtz and van Hove (HvH) theorem |27| . 
which is exact for all the interacting Fermi systems, independent of the type of the interaction 
between fermions. For the infinite NM, it is straightforward to see that the HvH theorem is 
satisfied on the HF level only when the in-medium NN interaction is density independent, 


i.e., when the RT is equal zero 


ft is, therefore, of high interest to assess the impact 


of the RT on the s/p potential in a mean-field study of the NM within the standard HF 
frame using a realistic density dependent NN interaction. Moreover, given the fact that the 
nuclear symmetry energy and nucleon effective mass are directly linked to the momentum- 
and density dependence of the single-nucleon potential [ 4 , 7|, it is highly desirable to have 
a method to take into account properly the density- and momentum dependence of the RT 
of the s/p potential on the HF level, which is the main motivation of the present study. It 
is further expected that, in a local density approximation, one should be able to include, in 
the same mean-field manner, the RT and momentum dependence of the density dependent 


NN interaction into the folding model calculation of 
which is being evaluated so far mostly on the HF level 


he nuc 


19 


20 


extension of the single-folding approach for the nucleon OP 


eon OP for the finite nuclei, 


22 


23]. Towards this goal, an 


19] has been done and applied 


to study the elastic neutron scattering on the lead target measured at the incident energies 
of 30.4 and 40 MeV 30]. 


II. EXTENDED HF FORMALISM FOR THE SINGLE-PARTICLE POTENTIAL 


In the present HF approach, we consider the homogeneous, spin-saturated NM of given 
neutron (p n ) and proton (p p ) densities, or equivalently, of given total nucleon density p = 
Pn + p P and neutron-proton (n/p) asymmetry 5 = (p n — p p )/p . Using the direct (t£*) and 
exchange (nf x ) parts of the (central) in-medium NN interaction v c , the total NM energy 
can be determined in the standard nonrelativistic HF scheme as E — E^ n + E pot , where the 
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kinetic and potential energies are 
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The s/p wave function \kor) is plane wave, and the summation in ([I])-© is done sepa¬ 
rately over the neutron (r = n) and proton (r = p) s/p indices. The nucleon momentum 
distribution n T {k) in the spin-saturated NM is a step function determined with the Fermi 
momentum k^ = (37T 2 p T ) 1 / 3 as 


n T [k) = 


1 if k ^ k!p 

0 otherwise. 


(3) 


m 


According to the Landau theory for the infinite Fermi systems 
in the NM is determined as 

r)F t> 2 k 2 

eAk) = = t T (k ) + U T {k) = — + U A k), 


26], the s/p energy e T (k) 


(4) 


which is the change of the NM energy caused by the removal or addition of a nucleon 
with the momentum k. The single-nucleon potential U T (k ) consists of both the HF and 
rearrangement terms 

UAk) = Ui UF \k) + Ui RT \k), (5) 

where U^ HF \k) — n T >(k')(kaT, fcW|u c |fc(7T, k'a'r')^ (6) 

k' ct't' 
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and U± RT \ k ) = - ^ n T1 (fci)n T2 (fc 2 


k\a\T\ k2CT2T2 

X ( feiUiTi, k 2 cr 2 r 2 
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kicnn, k 2 a 2 T 2 ) . 
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When the nucleon momentum approaches the Fermi momentum (k —> k ^), e T (kp' > ) deter¬ 
mined from Eqs. (HJ-(JZJ) is exactly the 
theorem 27}. Using the transformation 


ermi energy given by the Hugenholtz - van Hove 
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dn T (k ) 


dp T dkpd 


1 7T 2 
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where is the total volume of the NM in the momentum space, t 


le rearrangement term of 


the s/p potential U T at the Fermi momentum can be obtained 29j as 

n 


ur\k^kP)=i ' £ 


trr 2 ( 2jr ) 


TlT 2 
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k\T\.k 2 T 2 


dv r 


dk 
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n T1 (ki)n T2 (k 2 ) 


kiT 1 ,k 2 r 2 ) d 6 k l d 6 k 2 . 
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(9) 


At variance with the RT part, the HF part of the s/p potential can be readily evaluated at 
any momentum 

Rf F) jk) = 2 j~~~; [ n T '(k')(k,k'T'\v c \k, fe'r% d 3 k'. (10) 

^ (2vr) 3 J 

The spin components of plane waves in the HF and RT parts of the s/p potential are 
averaged out, and this results on the spin degeneracy factors 4 and 2 in the expressions (J9J) 
and (flop , respectively. We keep in mind that the s/p potential is determined consistently at 
each total NM density p. As a result, U T is a function of the total NM density p, the neutron- 
proton asymmetry <5, and the nucleon momentum k. For the spin-saturated NM, only the 
spin-independent terms of the central NN interaction are needed for the determination of the 


s/p potentials (l9lh(TTO. In 
(CDM3Y3 and CDM3Y6) 


he present work, we have used two density dependent versions 


18| of the M3Y interaction based on the G-matrix elements of 


the Paris NN potential in a oscillator basis 


12], Thus, the central part of the CDM3Yn 


interaction was used in the present HF calculation explicitly as 

v d ( ex) ( s ) _ F 0 (p)vJ 0 (EX) (s) + Fi(p)wJi (EX) (s)ti • t 2 , where s = |r*i * r 2 \. (11) 

The radial parts n^ EX ^(s) are kept unchanged as determined from the M3Y-Paris interaction 


12], in terms of three Yukawas 


D(EX), 
7 00 ( 01 ) ^ 


\^ v d(ex)^ ^exp(— R v s) 

r 00(0i) \ v ) R s 

U=1 


( 12 ) 


and the explicit Yukawa strengths and ranges can be found, e.g., in Table 1 of Ref. 20]. The 
density dependence of the interaction (fTT]) is assumed to have the same functional form as 


that introduced first in Ref. 


18j 


* 0 ( 1 ) (p) — Co(l)[l + CKo(l) ex p( — Po(i)P) + 7o(i)P]- 


(13) 
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FIG. 1. (Color online) Total NM energy per particle E/A at the different n/p asymmetries 5 
given by the HF calculation (P)-©, using the CDM3Y3 (lower panel) and CDM3Y6 (upper panel) 
interactions with their IV density dependence anew determined in the present work. The solid 
circles are the saturation densities of the NM at the different n/p asymmetries. 


The parameters of the isoscalar (IS) density dependence F 0 (p) were determined (l8f] to repro¬ 
duce the saturation properties of the symmetric NM and give the nuclear incompressibility 
K = 218 and 252 MeV with the CDM3Y3 and CDM3Y6 interactions, respectively. These 
interactions, especially the CDM3Y6 version, have been well tested in numerous folding 


model analyses of the elastic nucleon-nucleus 


19] and nucleus-nucleus scattering 21|, and 
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TABLE I. Parameters of the density dependence (1131) of the CDM3Yn interaction. The incom¬ 
pressibility K of the symmetric NM, the nuclear symmetry energy Sq and its slope L were obtained 
from the HF results (1301) at the saturation density po ~ 0.17 fm -3 . 


Interaction 

i 

Ci 

Ui 

Pi 

7 i 

K 

So 

L 





(fm 3 ) 

(fm 3 ) 

(MeV) 

(MeV) 

(MeV) 

CDM3Y3 

0 

0.2985 

3.4528 

2.6388 

-1.500 

218 

- 

- 


I 

0.2343 

7.6514 

9.7494 

6.6317 

- 

30.1 

49.6 

CDM3Y6 

0 

0.2658 

3.8033 

1.4099 

-4.000 

252 

- 

- 


1 

0.2313 

7.6800 

9.6498 

6.7202 

- 

30.1 

49.7 


the charge-exchange scattering to the isobar analog states 20, |23]. Like in Ref. 20], the 
parameters of the isovector (IV) density dependence Fi(p) were determined in the present 
work to reproduce the BHF results for the IV term of the microscopic nucleon OP in the 
asymmetric NM obtained by Jeukenne, Lejeune and Mahaux (JLM) (3l|, |32] - Because of the 
RT included into the extended HF calculation of the nucleon OP, the parameters obtained 


for the IV density dependence F\(p) are slightly different from those used earlier 2CJ, 123 ]. 


For convenience of the readers who are interested to use the CDM3Yn interaction in the HF 
or folding model calculation, the parameters of the density dependence are given explicitly 
in Tabled] 

The HF results for the total energy per particle E/A of the asymmetric NM are shown 
in Fig. CD One can see that the saturation density rapidly decreases with the increasing 
n/p asymmetry, and the pure neutron matter (<5 = 1) is unbound by the (in-medium) 
NN interaction. At the high NM densities, the E/A curves obtained with the CDM3Y6 
interaction are stiffer than those obtained with the CDM3Y3 interaction, and this is due 
to the higher nuclear incompressibility K given by the CDM3Y6 interaction. The behavior 
of the EOS of the asymmetric NM with the increasing n/p asymmetries shown in Fig. [I] 
is typical and similar to those observed earlier in the HF calculations of the NM using the 
different types of the in-medium (density dependent) NN interaction 15HI7]. 


Given the parametrization (ITT]) of the CDM3Y3 and CDM3Y6 interactions, the HF part 


























of the s/p potential can be explicitly obtained in terms of the IS and IV parts as 

Ul HF \p, 5, k) = F 0 ( P )ur Y \ P , k) ± F 1 (p)U? gs * ) {p, 5, k), (14) 

where U^ M3Y \p,k) = pj£ + J A 0 (r)vo*(r)j 0 (kr)d 3 r, 

and U[ M3Y \p,S,k) = pJ F 5 + J Ai(r)v F ^(r)j 0 (kr)d 3 r. 


Here A 0( i)(r) = p n ji(kPr) ± p P ji(k^’r), J^ {1) = j 
and ji(x) — 3ji(x)/x = 3(sina; — xcosx)/a: 3 . 




7 oo(oi) 


[r)d 3 r, 


(15) 


The (-) sign on the right-hand side of Eq. (TT4|) pertains to the single-proton (r = p) and 
(+) sign to the single-neutron (r = n) potentials. Because the original M3Y interaction 
is momentum independent, the momentum dependence of the HF potential ifTTl) is entirely 
determined by the exchange terms of Uf^ 3Y ^. 

Applying the HvH theorem, the RT of the s/p potential is also obtained explicitly in 
terms of the IS and IV parts, but at the Fermi momentum only 


v( rt )(p, 5 ,4 r) ) = ur \ P , 4 T) ) + ^i (rt) (p, 4 r) ) 


(16) 


where U^ RT \p, kp’) = 


(T ) _ 1 dFp (p) 


2 dp 


P 2 Jo + I Ai{r)v^{r)d 3 r 


EX/ 


and f/| RT \p, 5, kp 1 ) = 


(r)\ _ I^I(P) 


p 2 J?5 2 + / Ai(r)v»?(r)d 3 r 


EX/ 


2 dp 

One can see from Eq. (1T61) that the RT becomes zero if the density independent M3Y 


interaction is used in the HF calculation of t 
theorem is satisfied already on the HF level 


re^s/p potential. In such a case, the HvH 
At a given n/p asymmetry S, the IV 


part of the RT is the same for both the single-neutron (r = n) and single-proton (r = p) 
potentials, and that affects the simple representation of the nucleon OP by Lane |34l. l35| . 
where the IV parts of the neutron and proton OP are equal but with the opposite sign. 

In general, as seen from Eq. O, the RT of the s/p potential should be present at the 
arbitrary nucleon momenta. Microscopically, the momentum dependence of the RT was 
shown, in the BHF calculation of the NM [l, 3, 4|, to be due to the higher-order NN 
correlation, like the second-order diagram in the perturbative expansion of the mass operator 
or the contribution from the three-body forces etc. In the finite nuclei, the rearrangement 
effects in the nucleon removal reactions (which have about the same physics origin as the 
RT potential considered here) were shown 33j to be strongly dependent on the energy of 
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the stripping reaction, a clear indication of the momentum dependence of the RT potential. 
The question now is whether one can assess the momentum dependence of the RT of the s/p 
potential on the HF level. Making use of the factorized density dependence of the CDM3Y3 
and CDM3Y6 interactions, we suggest in the present work a rather simple method to include 
consistently a momentum dependent RT into the s/p potential in the same HF framework. 
An important constraint of this procedure is that adding a realistic momentum dependent 


RT to the HF potential should improve the agreement o 


the calculated nucleon OP in the 


NM with the empirical data. It has been shown jl3, 24] that the momentum dependence 
of the HF potential (fTTD could account fairly well for the observed energy dependence of 
the nucleon OP after a slight adjustment of the interaction strength at the high energies. 
Therefore, we adopt phenomenologically a momentum dependent RT of the s/p potential in 
the functional form similar to (1T4|) as 

U^(p,6,k) = AF 0 (p)ul i mY) (p,k) + AF l (t>,6)ui M3Y \f,,6,k), (17) 


where the (momentum-independent) rearrangement contributions to the IS and IV den¬ 
sity dependences of the CDM3Yn interactions are determined consistently from the exact 
expression (TT6li of the RT at the Fermi momentum as 

ur\p,kP) U^MkP) 


A F 0 (p) = 




Ob 


and AFi(p, 5) = 


ui M3Y \ P ,s,k^k ( ;>) 


Oh 


(18) 


Consequently, the total s/p potential is determined in the present HF approach as 


Ur(p,6, k) = Uo(p, k) ± Ui(p, 5, k ) 

= [F 0 (p) + AF 0 (p))£/< M3Y) (Ai) 

± (FiW ± AF 1 (p,*)]E7 1 (M3Y) (M*0, (19) 

where (-) sign pertains to t = p and (+) sign to r — n. Thus, the momentum dependence 
of the total s/p potential U T is determined by that of the exchange terms of Due to 

the presence of the RT, the absolute strength of the IV term of the single-proton potential 
is not equal that of the single-neutron potential. One can see from the expressions (fT71h(fT9l) 
that the rearrangement effects actually result on a modification of the IS and IV density 
dependence of the central interaction (fill) . v c —* v c + Av c , so that the total s/p potential 
can be estimated in the standard HF scheme as 


U T (p, S, k ) = n T >(k')(kar, k'a'r'\v c + Au c |fccrr, k'a'r')^. (20) 

k'cr'r 1 


10 







Nucleon OP in the NM and the s/p potential at high momenta 


In the NM limit, the nucleon OP is determined as the (mean-held) interaction potential 


between the nuc 
filled Fermi sea 


eon incident on the NM at a given energy E and the bound nucleons in 


24]. In general, the nucleon OP contains both the IS and IV parts [34, 


he 


351 


like the total s/p potential (fl9j) . Given a strong dominance of the IS term of the nucleon 


OP 


19 


22j, one needs first to explore the IS term of the nucleon OP predicted by t 


he HF 


calculation of the NM. Applying a continuous choice for the nucleon s/p potential 


36 ] at 


the positive energies E, we obtain in the HF scheme the nucleon OP in the symmetric NM 
jl3, 24] as 


Uo(p, E) = Ul™\p, E) = F 0 (p)p 


Jf + I ]Jkrr>J a (kiH,p,)r)if^{r)d'r 


( 21 ) 


Here k(E,p) is the (energy dependent) momentum of the incident nucleon propagating in 
the mean held of the nucleons bound in the NM, and is determined as 


/ 2777 / 

k(E,p) = — [E - U 0 (p,E)}, with E > 0. 


( 22 ) 


It is easy to see that k(E , p) > fcp and t/jg F is just the high momentum part of the isoscalar 
term of the HF potential (1I4|) . Based on the above discussion, the total nucleon OP in the 
NM should also have a contribution from the RT added 


Uo(p, E) = uiP\p, E) + u‘ s RT> (p, E) 

= PoW + A F„(p)]t/' M 3 Y) (M(£,p)), 


(23) 


where the density dependence AF 0 (p) of the RT is determined by the relation (US]) . 

The total nucleon OP (]23l) evaluated at the saturation density po of the symmetric NM 


using the CDM3Y6 interaction are compared with the empirical data 


3 - 39 ] 


in Fig. [2] Al¬ 


though the inclusion of the RT significantly improved the agreement of the calculated Uq 
with the data at the lowest energies, it remains somewhat more attractive at the higher 
energies in comparison with the empirical trend. Such an effect is easily understood in light 


of the microscopic BHF results for the nucleon OP 


36], where the energy dependence was 


shown to come not only from the exchange part, but also from the direct part of the micro¬ 
scopic OP because of the energy dependence of the Brueckner G-matrix. That is the reason 
why a slight linear energy dependence has been introduced into the CDM3Y6 interaction 


11 

























FIG. 2. (Color online) Energy dependence of the nucleon OP in the symmetric NM (evaluated at the 
saturation density po with and without the RT u sing the CDM3Y6 interaction) in comparison with 


the empirical data taken from Refs. m (circles), 


• 137] (cii 


38j (squares) and 


311] (triangles). The momentum 


dependent factor g(k ) has been iteratively adjusted to the best agreement of the calculated nucleon 
OP (124|) with the empirical data (solid line). 


18l. |l9j, in terms of the g(E) factor. To be consistent with the momentum dependence of 
the s/p potential under study, instead of the g(E ) factor, we scale in the present work the 
CDM3Yn interaction (ITT]) by a momentum dependent function g(k(E, p)), and iteratively 
adjust its strength to the best agreement of the (HF+RT) nucleon OP obtained at the 
saturation density p 0 with the empirical data, as shown in Fig. [2] As a result, 

(M3Y) 


Uo(p, E ) = g(k(E, p)) (F„(p) + AF„(p)](7'“ v ’ (p, k(E, p)), 


(24) 

where k(E , p ) is determined self-consistently from Uq(E, p ) by Eq. (J22J). At variance with the 
g{E) factor fixed by the incident energy 18|, ll9|], g{k{E 1 p)) is now a momentum dependent 
function (see Fig. [3]), carrying the important signature of the momentum dependence of the 
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FIG. 3. Momentum dependent scaling factor g(k) obtained with both the CDM3Y3 and CDM3Y6 
interactions from the best (HF+RT) fit of the nucleon OP (1241) to the empirical energy dependence 


39]. 


of the nucleon OP l39f]. The points are the numerical results that are well reproduced by a cubic 
polynomial (solid line). 


nucleon mean-field potential. Numerically, the obtained g(k ) function is nearly identical 
for both the CDM3Y3 and CDM3Y6 interactions, and it can be considered as the explicit 
momentum dependence of the CDM3Yn interaction that allows the incident nucleon to feel 
the nucleon mean-field potential during its interaction with the nucleons bound in the NM. In 
this sense, such a momentum dependence is of a similar nature as the momentum dependence 
of the G-matrix in the microscopic BHF study of NM, which is determined self-consistently 
through the momentum dependence of the s/p energies embedded in the denominator of the 
Bethe-Goldstone equation 2, ^]. The technical difference here is that the /^-dependence of 
g(k) has been determined from the best fit of the calculated s/p potential (I23|) at positive 
energies with the observed energy dependence of the nucleon OP. It can be seen in Fig. [3] 
that g{k) becomes smaller unity at k > 1.6 fm _1 only. Consequently, the obtained g(k) 
function is used further in the extended HF calculation to adjust the high-momentum part 
of the (HF+RT) s/p potential. 
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FIG. 4. Density dependence of the IV part of the neutron OP in the pure neutron matter at 
the energy E = 1 MeV, evaluated with and without the RT. The parameters (1131) of the IV 
density dependence F\(p) of the CDM3Y6 interaction have been iteratively adjusted to give the 
best agreement of the HF+RT result with that of the BHF calculation by the JLM group (solid 


circles) 


31, 3 - 


Situation is quite different concerning the IV term of the nucleon OP because there are 
no systematic (energy dependent) empirical data available, like those discussed above for 
the IS potential. However, it is well established from numerous optical model analyses of 
the elastic nucleon scattering that the absolute strength of the IV term of the nucleon OP is 


much weaker than that of the IS 




38]. and 


dominantly determined by that of the IS term 


he energy dependence of the nucleon OP is 
. Consequently, the momentum dependent 


function g(k) determined above for the IS part should be a reasonable approximation for 
the IV part of the nucleon OP in the HF+RT calculation 


V\(p, 5, E) = g(k(E, 6, p)) [fi(p) + Af\(p, 5)]C< M3Y| (p, S, k(E, S,p)), 


(25) 


where k(E , S, p) is determined self-consistently from the total nucleon OP in the asymmetric 
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NM as 


/ 2 rn 

k(E,i, f ,)=^j—\E-U„(f,,E)^U l (p,£)]. 


( 26 ) 


Thus, the NN interaction in the t.t channel of the central force (ITTj) is also influenced by 
the nucleon mean-field potential through the momentum dependent g(k ) function. Like in 


Refs. 


20 


23], parameters of the IV density dependence Fj (p) of the M3Y-Paris interaction 


have been iteratively adjusted in the HF+RT calculation to achieve a good agreement of 
the IV potential (1251) with the microscopic nucleon OP in the asymmetric NM given by the 
BHF calculation done by the JLM group [h, 32], with the contribution of the RT properly 
included 40]. At variance with the previous studies where F\(p) has been determined 


separately at each considered energy 


20 


set of the parameters of F\ (p) which were determined from the 
to the corresponding JLM results at the lowest energies 


231 ], we have used in the present work a unique 

3 est fit of the IV potential 


32 ]. where g(k) ~ 1 (see 


Fig. [4]). Together with the IS density dependence F 0 (p) determined earlier jl8|, the newly 
determined IV density dependence Fi(p) of the CDM3Yn interactions have been used in 
the present work to calculate the total NM energy and the single-nucleon potential. The 
absolute strength of Fi(p) in the extended HF calculation was scaled by a factor of 1.3, 
which was found necessary in the folding model analysis of the ( p , n ) scattering to the isobar 
analog states |20]. As a result, the nuclear symmetry energy S'(po) at the saturation density 
given by the present HF calculation is very close to the empirical value of about 30 MeV 
(see also next Section). The final values of the parameters of Fo(p) and Fi(p) are given with 
the most important NM properties in Table UJ 

The total (density- and momentum dependent) s/p potential is now determined as 


U r {p, <5, k) = g(k)[U 0 (p, k ) ± £A(p, <5, k)], (27) 

where U 0 (p, k ) and Ui(p, 5, k ) are determined by using Eq. (fT9jl and the same function g{k) 
as that shown in Fig. [3] The total s/p potentials obtained at the saturation density p 0 ~ 0.17 
fm -3 from the extended HF calculation (l20jl of the symmetric NM and the pure neutron 
matter using the CDM3Y6 interaction are shown in Figs. [5] and |6] respectively. One can see 
that the RT is largest at the nucleon momenta close to zero, which correspond to nucleons 
deeply bound in the NM (k <C kp). The RT steadily decreases with the increasing nucleon 
momentum, and the decrease of the RT becomes faster when the high-momentum part of the 
s/p potential is scaled by the g(k) function determined [see Eqs. (T24j) and (l25j) ] to reproduce 
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FIG. 5. Momentum dependence of the total s/p potential in the symmetric NM at the saturation 
density po, with the explicit contributions from the RT and HF parts. The upper panel shows 
the results of the extended HF calculation (1231) . and the lower panel shows the total s/p potential 
(1241) . with the high-momentum part corrected by the g(k) function determined from the observed 
energy dependence of the nucleon OP. 


the observed energy dependence of the nucleon OP. In this case, the s/p potential reaches 
zero and changes sign at the momentum k « 3.3 fm -1 that corresponds to the nucleon OP 
at the incident energy around 220 MeV, where the empirical (Schrodinger-equivalent) U op 
was shown to become repulsive [39]. Such a momentum dependence of the s/p potential 


agrees well with that predicted by the microscopic BHF calculation (see, e.g., Figs. 6 and 
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FIG. 6. The same as Fig. [5] but for the total s/p potential in the pure neutron matter. 


7 of Ref. [3j). A decrease of the RT with the increasing nucleon momentum (or energy) 
also agrees with the observed energy dependence of the rearrangement effect in the nucleon 
removal reactions 33|. In the relative strength, the RT contributes to about 20 ~ 30 % of 


the total strength of the s/p potential at p = p 0 over a wide range of the nucleon momentum 
and is, therefore, a very clear manifestation of the HvH theorem [ 27 I ] . 


The dependence of the neutron and proton s/p potentials on the n/p asymmetry (with 
the explicit contributions from the RT and HF parts) is shown in Fig. [71 and one can see 
that the RT is the same for both the neutron and proton s/p potentials. In the present 
extended HF scheme, such an equality is exactly obtained from the RT given by the HvH 
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FIG. 7. (Color online) Contributions of the RT and HF parts to the total neutron and proton 
s/p potentials evaluated at the different n/p asymmetries 5 and nucleon momenta k = kp \ using 
the CDM3Y6 interaction. The results obtained at the saturation density p = po are shown in the 
upper panel, and those obtained at the density p = 2po are shown in the lower panel. 


theorem (fTRlb At a given NM density, the RT decreases slightly with the increasing n/p 
asymmetry 5. As can bee seen from the lower panel of Fig. [TJ the repulsive contribution of 
the RT becomes much stronger at the high density p = 2p 0 , with the relative strength up to 


70 % of the HF term. Such a behavior of the RT is well expected, given the 


correlations and the three-body forces as the physics origin of the RT [1, 3|, which become 


righer-order NN 
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much more substantial with the increasing NM density. The results of our extended HF 
calculation shown in Fig. [7] also agree well with those of the recent BHF calculation of the 
asymmetric NM by Vidana 4l|, using the Argon NN interaction (Vi 8 version 42] with the 
Urbana three-body force). 

To conclude this section, a simple and consistent method has been developed to account 
effectively for the momentum dependence of the RT of the s/p potential in an extended 
HF calculation using the CDM3Yn density dependent interactions (fill) , based on the exact 
expression of the RT given by the HvH theorem at each NM density and the empirical energy 
dependence of the nucleon OP observed over a wide range of energies. 


III. SINGLE-NUCLEON POTENTIAL AND THE SYMMETRY ENERGY 

Given the importance of the nuclear symmetry energy for the nuclear astrophysics studies, 
especially, its vital role in the determination of the EOS of the asymmetric NM [§|, 9], we 
focus in the present section on the connection of the s/p potential in the NM with the nuclear 
symmetry energy, which has been widely discussed in Refs. 4, 7]. The nuclear symmetry 
energy S(p) is normally defined with an expansion of the total NM energy per particle over 
the n/p asymmetry <5 as 


j(p,S) = j(p, i = o) + S(p)S 2 + 0(6‘) + 


(28) 


For S < 1, the contribution of 0(<5 4 ) and the higher-order terms in Eq. (1281) was proven to 


be quite small j3|, [l5| and can be neglected (the so-called parabolic approximation). For the 
pure neutron matter (5 = 1) the convergence of the series (|28l) is slower, and the parabolic 


law becomes less accurate. Therefore, a more general definition of the nuclear symmetry 
energy as the energy required per particle to change the symmetric NM into the pure neutron 
matter is also used in the mean-held studies 


S(P) = 3 = 1)- ^(p, 5 = 0). 


The nuclear symmetry energy S(p) can also be expanded [44] around p 0 as 


sm = s 0 + %(?—^)+^(?—^'i +.. 

3 V Po J 18 V Po 


(29) 


(30) 


where L and K sym are the slope and curvature parameters of the symmetry energy at po- 


While the curvature parameter K sym is still poorly known, it has been shown recently by Li 
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FIG. 8. (Color online) Nuclear symmetry energy S(p) at the n/p asymmetries 5 = 0.6 and 1 
obtained with the CDM3Y6 interaction, using the parabolic approximation (1281) and general defi¬ 
nition (1291) . The shaded (magenta) region marks the empirical boundaries implied by the analysis 
of the isospin diffusion data and the double ratio of the neutron and proton spectra observed in 
the HI collisions [43l . 44]. The square is the empirical S value implied by the structure study of the 


GDR 


451 ]. and the triangle is that established from the analysis of the terrestrial nuclear physics 


experiments and astrophysica] 
calculations by Akmal et al 


observations 


[ 47 ]. 


The circles and crosses are results of the ab-initio 


48] and Gandolfi et al. 


491 ]. respectively. 


and Han j47j that quite a robust constraint for both Sq and L values has been established 
based on several tens analyses of the terrestrial nuclear physics experiments and astrophysical 
observations, which give Sq ~ 31.6±2.7 MeV and L & 58.9±16.0 MeV. With the parameters 
of the IV density dependence of the CDM3Yn interaction anew determined in the present 
work, the values of So ~ 30.1 MeV and L « 49.7 MeV given by the HF calculation are well 
within this empirical range. 


Applying the HvH theorem to calculate^t 
responding Fermi momenta, one obtains 


re neutron and proton energies (j4j) at the cor- 
7] the nuclear symmetry energy directly from 
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the difference between the neutron and proton Fermi energies as 


tn(4 n) ) - tp(k^) + u n (p , 5, 4 n) ) - U p {p, S , 4 P) ) = 4 S(p)S + 0(S 3 ) + ... (31) 


Thus, Eq. (T5TT) is simply the parabolic law in the s/p energy representation. As discussed in 
Sect. M at given p and 5 values the contribution of the RT to the s/p potential (TT6T) is the 
same for both the neutron and proton s/p potentials. Therefore, the RT contribution to the 
nuclear symmetry energy S(p) through the difference between the neutron and proton s/p 
potentials in Eq. (l3Tj) is canceled out. The same conclusion was also drawn in the recent 
BHF study of asymmetric NM by Vidana 4l|. Using the single-nucleon potentials given 
by the present HF calculation, we obtained exactly the same S(p) from both Eq. (j3TT) and 
the expansion (128]) . neglecting the higher-order terms 0(5 3 ) and 0(<5 4 ), respectively. These 
results are compared with those given by the general definition (1291) in Fig. [HJ One can see 
that the parabolic approximation is reasonable only for 6 < 1, and it becomes poorer for 
pure neutron matter. As a result, the higher-order terms on the right-hand side of Eqs. (|28li 
and (T3TT) need to be taken into account for very neutron-rich matter, as done recently by 
Chen et al. jf|. 

The main method to probe S(p) obtained with a chosen in-medium NN interaction is to 
probe this interaction in the analysis of heavy-ion (HI) collisions 43, |44j or in the structure 


studies of nuclei with large neutron excess |45l. 1461 ] . Based on the constraints implied by 
such studies, extrapolation is often made to study the low- and high-density behavior of 
nuclear symmetry energy. For the illustration, we have compared in Fig. E] the HF results 
given by the CDM3Y6 interaction for S(p) with the empirical data 143H47] and the results 


of 


he ab-initio calculations of the asymmetric NM by Akrnal et al. [48j] and Gandolh et al. 
49]. Around the saturation density p 0 the symmetry energy S7 given by the HF calculation 


is in a very good agreement with the empirical value 46j, (47]. In the low-density region 
(p ~ 0.3 ~ 0.6p 0 ) the empirical boundaries for S(p) deduced from the analysis of the isospin 


diffusion data and double ratio of neutron and proton spectra data of HI collisions [43l , 144 ] 
do enclose our HF result. At the density p ~ 0.1 fm -3 , the HF result also agrees well with 
the empirical value deduced from the structure study of the Giant Dipole Resonance (GDR) 
in heavy nuclei 45]. At supranuclear densities, where the reliable empirical data are still 


absent, our HF results follow closely those of the ab-initio calculations 


48 


49]. 


Although the explicit contribution of the RT to the nuclear symmetry energy, through 
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FIG. 9. (Color online) Nuclear symmetry energy S(p) given by the HF calculation using the new 
parameters of the IV density dependence F\ (p) of the CDM3Y6 (upper panel) and CDM3Y3 (lower 
panel) interactions. These same HF results obtained with the old parameters of F\(p) determined 


in Ref. 


201 ] are shown as the dashed lines. Other symbols are the same as in Fig. [8] 


the difference between the neutron and proton s/p potentials f[3T]) . is canceled out. The vital 
role of the rearrangement effects in the HF calculation of the nuclear symmetry energy is 
well illustrated in Fig. [9l where S(p) values obtained from the HF calculation using the new 
parameters of the IV density dependence F\ (p) of the CDM3Yn interactions are compared 
with those obtained with the old parameters of F\ (p ), determined in Ref. 20] by fitting the 
IV part of the (HF only) s/p potential to the IV part of the microscopic JLM potential. One 
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FIG. 10. (Color online) Nuclear symmetry energy S(p) obtained at the n/p asymmetry 5 < 1 
from the (HF+RT) s/p energies with the CDM3Y6 (upper panel) and CDM3Y3 (lower panel) 
interactions, using the parabolic approximation (1311) . The solid (dashed) curves show the results 
obtained with (without) the modification of the high-momentum part of the s/p potential by the 
g(k) function determined from the observed energy dependence of the nucleon OP. Other symbols 
are the same as in Fig. [8] 

can see in Fig. [2] that S(p) given by the old IV density dependence agrees with the empirical 
data and ab-initio results at NM densities up to about 2/?o only. Using the new parameters 
of Fi(p) obtained in the present work by fitting the IV part of the (HF+RT) s/p potential 
to the JLM potential (see Fig. QJ, the calculated S(p) values now follow closely those given 
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by the ab-initio calculations over a wider range of the NM densities, up to p ~ 4po- 

With the Fermi momentum Jzf becoming larger at high NM densities, a realistic momen¬ 
tum dependence of the s/p potentials should be helpful in constraining the nuclear symmetry 
energy S(p) at the supranuclear densities, based on the relation (I3ip . The symmetry en¬ 
ergies S(p) given by the difference between the neutron- and proton s/p energies (T3TT) at 
the n/p asymmetries 6 < 1 are shown in Fig. [TUI with and without the modification of the 
high-momentum part of the (HF+RT) nucleon OP by the g(kp ) function determined from 
the observed energy dependence of the nucleon OP. One can see that the modification of the 
high-momentum tail of the (HF+RT) s/p potential results on a slightly softer slope of the 
nuclear symmetry energy at the high NM densities that leads, in turn, to a good agreement 
of the HF results with those of the ab-initio calculations over a much wider range of the NM 
density, up to p ~ 5 ~ 6/+- This result is, thus, complementary to the recent efforts by Li 


et al. 


50] to determine the nuclear symmetry from the optical model analysis of the elastic 


neutron-nucleus scattering over a wide range of energies. 


IV. NEUTRON-PROTON EFFECTIVE MASS SPLITTING 


As discussed in Sec. m due to the finite range of the Yukawa functions in the radial 
part of the CDM3Yn interaction (TTTTl . the s/p potential depends explicitly on the nucleon 
momentum k through its exchange term that implies a nonlocal single-nucleon potential in 
the coordinate space. At the high nucleon momenta, the momentum dependence of the s/p 
potential is further modified by the g{k) function implied by the observed energy dependence 
of the nucleon OP. An important quantity associated with the momentum dependence of 
the nucleon s/p potential is the nucleon effective mass m*, defined within the nonrelativistic 
mean-held formalism as 


m*{p,5) 


m 


1 + 


m dU T (p,S,k) 


h 2 kP 9k 


S T ) 


-1 


(32) 


where m is the free nucleon mass. Thus, the nucleon effective mass m* describes the nonlo¬ 
cality of the mean-held potential felt by a nucleon propagating through the nuclear medium. 
As such, the nucleon effective mass is directly linked to many nuclear physics phenomena, 
like the dynamics of HI collisions, damping of giant resonances, temperature prohle of the hot 


stellar objects and neutrino emission therefrom 


5l|. For the asymmetric NM, the relative 
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FIG. 11. (Color online) Density dependence of the neutron effective mass (1321) at the n/p asym¬ 
metry 5 = 0, 0.3 and 0.6, obtained from the s/p potential given by the CDM3Y6 interaction. The 
HF and HF+RT results are shown in the upper and lower panels, respectively. The solid circles 
are the effective-mass values determined at the saturation density po- 


difference between the neutron and proton effective masses 


m* n _ p {p,6) = 


m 


(33) 


is widely discussed [d?, 521 as the neutron-proton effective mass splitting , which is closely 
related to the nuclear symmetry energy S(p) and its slope parameter L j^]. The n/p 
effective mass splitting was also suggested to affect the neutron/proton ratio during the 
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FIG. 12. (Color online) The same as Fig. [II] but for the proton effective mass. 


stellar evolution, and the cooling of protoneutron stars etc. 52]. 


In the present work we do not intend to explore this issue over such a wide scope, but 
focus briefly on the effect of the RT in the determination of the nucleon effective mass 
within the extended HF formalism as well as the effect caused by the modification of the 
high-momentum part of the s/p potential to match the observed energy dependence of the 
nucleon OP. The obtained neutron and proton effective masses are shown in Figs, fill and 
Figs. H21 respectively, and one can see that the RT enhances m* substantially at high NM 
densities. A similar behavior of the nucleon effective mass can also be seen in the results of 
a recent microscopic BHF calculation by Baldo et al. [51], where the RT originating from 
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FIG. 13. (Color online) The same (HF+RT) results for the neutron effective mass as in the lower 
panel of Fig. [TTJ obtained with (lower panel) or without (upper panel) the modification of the 
high-momentum part of the s/p potential by the g(k) function implied by the observed energy 
dependence of the nucleon OP. 


three-body force drastically enhances the nucleon effective mass at the high NM densities 
(see, e.g., Fig. 1 of Ref. |5l|). Like found in Sec. IIIII for the nuclear symmetry energy, the 
modification of the high-momentum part of the (HF+RT) s/p potential by the g{k) function 
implied by the observed energy dependence of the nucleon OP changes slightly the slope of 
the density dependence of the nucleon effective mass at the high NM densities (see Figs. [13] 
and fill) . Although the nucleon effective mass is still poorly known at the high NM densities 
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FIG. 14. The same as Fig. [13] but for the proton effective mass. 


and/or the large n/p asymmetries, the empirical m* jm value in the symmetric NM is known 


to be about 0.73 [53j at the saturation density p 0 . In a good agreement with that empirical 
data, our HF and HF+RT results obtained with the CDM3Y6 interaction (see Fig. [T3|) give 
m*/m ~ 0.734 and 0.755, respectively. The modification of the high-momentum part of 
the (HF+RT) single-nucleon potential by the g[k) function reduces this result slightly to 
m*/m ~ 0.737 at p = p 0 . 


With quite a good agreement of the calculated nucleon effective mass with the empirical 
data at 5 = 0, it is of interest to consider further the 5-dependence of the m*/m value. 
The HF results obtained at the saturation density p 0 are shown in Fig. [151 Given the 
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FIG. 15. (Color online) Neutron and proton effective masses (upper panel) and their splitting 
(lower panel) obtained at p = po and the different n/p asymmetries 5. Both the HF and HF+RT 
results show a well defined linear 5-dependence of the n/p effective mass splitting (1331) . 


equal RT contribution to each of the neutron- and proton s/p potential, the derivative with 
respect to the nuc;eon momentum gives about the same RT contribution of around 10% 
to the neutron- and proton effective masses (upper panel of Fig. [T5]) . and the n/p effective 
mass splitting (lower panel of Fig. H5]) is rather weakly affected by the rearrangement effect. 
On can see that both the standard HF and extended HF+RT results show a well defined 
linear 5-dependence of the n/p effective mass splitting (J33J). Our final HF+RT results 
give m*_ p (p 0 ,5) ~ (0.26 ± 0.01)5 which is well inside the empirical boundary m*_ p (p 0 , 5) ~ 
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FIG. 16. 5-dependence of the neutron-proton effective mass splitting obtained at p = 2/?o- The 
HF+RT results were obtained with (solid line) or without (dashed line) the modification of the 
high-momentum part of the s/p potential by the g(k) function determined from the observed energy 
dependence of the nucleon OP. 


(0.27+0.35)5, established from the analyses of the terrestrial nuclear physics experiments and 


astrophysical observations |47j. Our result is, however, lower than that estimated recently, 
ml (pn, 5) ~ (0.41 ± 0.15)5, from the phenomenological (isospin-dependent) nucleon OP 
52t ] determined from the extensive optical model analysis of a large data set of the nucleon 


elastic scattering. The uncertainty of 0.015 in the HF+RT result is not statistical error 
associated with the uncertainties of the model ingredients, but the uncertainty in adopting 
a linear 5- dependence of the n/p effective mass splitting (f+Tll in Fig. [T5l This fact indicates 
simply that at 5 < 1, when the parabolic approximation (1281) and (l3Tf is reasonable for 
the nuclear symmetry energy, the first-order symmetry term in the expansion of the single¬ 
nucleon potential over 5 (see, e.g., Ref. B) contributes overwhelmingly to the determination 
of the m*„ p (p, 5) value. 

With larger Fermi momentum Jzf at the high NM densities, the modification of the 
momentum dependence of the s/p potentials by the g(kp) function determined from the 
observed energy dependence of the nucleon OP should be taken into account in the HF+RT 
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calculation of the n/p effective mass splitting. One can see in Fig. U6lthat 5) becomes 

larger at p = 2po, and is enhanced further by more than 50% when the modification of the 
high-momentum part of the s/p potential is taken into account. The behavior of the nucleon 
effective mass (|52j) and the n/p effective mass splitting (fool) at the supranuclear densities 
is still poorly known, and the difference in the HF+RT results for the n/p effective mass 
splitting, m* n _ p (p, 5) ~ (0.48 ± 0.01)5 at p = 2 p 0 compared to (0.26 ± 0.01)5 at p = po, is 
quite significant and should be of interest for the nuclear astrophysical studies. 


V. FOLDING MODEL OF THE NUCLEON OPTICAL POTENTIAL 


Given the substantial rearrangement effects to the nucleon OP found above in the ex¬ 
tended HF calculation of the NM, it is of high interest to study these effects in the many-body 
calculation of the nucleon OP of the finite nuclei. As such, the folding model (more precisely 
biding model) has been proven to be a very effective tool to estimate the nucleon 


OP 


sine 

le-f 

19, 

54, 


55j . It can be seen from the basic folding formulas that this model generates the 


irst-order term of the microscopic OP defined in FeshbaelTs formalism of nuclear reactions 


56], based on the nucleon degrees of freedom. The success of the single-folding approach 


in the description of the elastic nucleon-nucleus scattering data measured for the targets in 
the different mass regions suggests that the first-order term of the FeshbaelTs microscopic 
OP is indeed the dominant part of the nucleon OP. 

In the single-folding approach, the central OP of the r-kind nucleon incident on the target 
A at the energy E is evaluated as a HF-type potential [19], using an appropriate, energy- 
and density dependent effective NN interaction u c (p, E) 


U‘° li (E,R) = X)[(fc T ,jk D (p, E)\k T ,j) + (k T ,j\vf x (.p,E)\j,k T )}, 

jeA 

= j\v c {p, E)\k T , j) A , 

jeA 


(34) 


where | j) is the single-nucleon wave function of the j —th nucleon of the target. The anti- 
symmetrization A of the nucleon-nucleus system is done by taking into account explicitly 
the knock-on exchange effects. As a result, the exchange term of U T becomes nonlocal in 
the coordinate space ( 55 )]. An accurate local approximation is usually made by treating 


the relative motion locally as a plane wave 


19 


54] . and the local energy dependent folded 
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potential (j34|) is obtained as an explicit function of the nucleon-nucleus distance R and local 
momentum k(R) of the incident nucleon. 


At the low incident energies, the pair-wise interaction between the incident nucleon and 
the nucleons bound in the target can induce certain rearrangement of the s/p configurations 
of the target nucleons that can be observed experimentally in the nucleon removal reactions 


33j. In terms of the nucleon-nucleus interaction, such a rearrangement effect is expected 


to affect also the shape and strength of the nucleon-nucleus OP (1341) . constructed in the 
folding model on the HF level from the s/p wave functions of the target nucleons. On the 
other hand, if we make a local density approximation (LDA) for the s/p potential (12UD by 
replacing the plane waves | k'a'r') by the s/p wave functions | j) of the target nucleons, then 
the resulted potential is just the HF-type folded potential (l34j) added by a rearrangement 
term, through the contribution of the Av c term. Thus, the total (central) nucleon OP now 
becomes 


Ur(E,R) = U? F (k,R) + U? T (k,R) = J2(k T J\v c (p,E) + Av c (p)\k T ,j) A , (35) 

j&A 


with I/ RT originated from the explicit density dependence of the CDM3Yn interaction (JTTJ) . 
In the standard HF formalism for the nucleon mean-held potential, the HF potential for an 

n 

unbound (scattering) nucleon has been assumed years ago [57], |58j as the first-order term of 
the microscopic OP for the low-energy elastic nucleon scattering, in about the same scheme 


as that of the Feshbach’s microscopic OP 


561 ]. Applying the variational HF method to 


obtain the scattering equation for the nucleon-nucleus system, a rearrangement contribution 
to the HF potential appears naturally if the effective NN interaction is density dependent 


591 ]. Given the rearrangement contribution to the density dependence of the CDM3Yn 


interaction determined above in the HF study of the NM, it is now possible to consistently 
include the RT into the folding calculation of the nucleon-nucleus potential 


Making explicit the neutron and proton s/p wave functions in Eq. (1351) and using the same 
isospin-dependent CDM3Yn interactions (ITT]) . we obtain the real nucleon OP explicitly in 
terms of the (energy dependent) isoscalar and isovector parts as 


Ur(E, R) = U m (E , R) ± U W (E , R), 


(36) 
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FIG. 17. (Color online) Real n+ 208 Pb optical potential at the neutron incident energies of 30.4 
and 40 MeV predicted by the single-folding calculation (j37|) - ([38]) using the CDM3Y6 interaction, 
with and without the contributions of the RT and the scaling g(k(E,Rj) function. CH89 is the 
phenomenological OP taken from the global systematics by Varner et al. 
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where (-) sign pertains to the proton OP and (+) sign to the neutron OP. 


U m (E,R) = g(k(E,R)) j [F„(p(r)) + AF„(p(r))J {Mr) + ftMWM 
+ [p„(R,r) + p r (R,r)]v$f(s)j 0 (k(E, R)s)}d?r, 

U„{E,R) = g(k(E,R)) J [fi(p(r» + AFi(p(r),T,(r))] {[ftM - ft W]' 
+ {Pn(R, r) - p p (R, r)]v$f(s)j 0 (k(E, R)a) J<i 3 r. 




(37) 

(38) 
and 


Here, p T (r,r') is the nonlocal s/p density matrix of the target with p T (r) = p T (r,r ) 
k(E, R) is the local momentum of the incident nucleon determined self-consistently from the 
total real nucleon OP as 


k(E,R) = J^[Ep. m .-U r (E,R)]. 


(39) 


Similarly to the expression of the nucleon OP in the NM limit (1271) . the folded nucleon- 
nucleus potential (1361) also depends explicitly on the momentum (l3fll) of the incident nucleon 
through the ( localized ) exchange term. To estimate properly the contribution of the re¬ 
arrangement effects, the IS and IV density dependences of the RT in Eqs. (I3?j) -([38 j) are 
determined in the LDA from the A F 0 (p) and AF\(p, <5) values given by the HF study of the 
NM (JTSjh accurately interpolated for the local nuclear density p = p[r) = p n (r) + p p (r) and 
the neutron-proton asymmetry 5 = T z (r ) = \p n {x) — Pp( r )\/p( r )- 

It can be seen from Eqs. (I37|l - (l38|) that the energy dependence of the real nucleon OP (l36|) 
is determined entirely by the local momentum of the incident nucleon k(E,R) appearing 
in the exchange potential as well as in the g{k{E 1 R)) function. Given the g(k ) function 
determined above in the HF study of the NM based on the observed energy dependence of 
the nucleon OP, each local scaling factor g(k(E, R)) of the folded potential is interpolated 
from the g{k) function (see Fig. [3]) at the local momentum k = k(E,R). As a result, 
g(h(E, R)) can now be considered as the explicit, energy (or momentum) dependence of the 
density dependent CDM3Yn interaction (fill) , locally consistent with the nucleon mean-held 
potential (1361) . This is an essential improvement of the present formulation of the single¬ 
folding model, compared to the earlier applications of the folding model (see, e.g., Ref. jl9|) 
where a constant factor g(E) « 1 — 0.0026F was used to scale the CDM3Y6 interaction. 
Because of the self-consistent determination of the g(k(E,R)) function and contribution of 
the RT through AF 0 (p(r)) and A Fi(p(r),T z (r)), the single-folding calculation (137|1 - 
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FIG. 18. (Color online) Elastic n+ 208 Pb scattering cross sections at the neutron incident energies 
of 30.4 and 40 MeV obtained_ with the real OP’s shown in Fig. [TTJ in comparison with the data 
measured by DeVito et al. 


global systematics CH89 


301 ]. The imaginary and spin-orbital parts of the OP taken from the 


381] were used in the OM calculation. 


becomes more cumbersome and time consuming compared with the earlier version jl9|] of 
the folding model. 

Although a comprehensive folding model study of the elastic nucleon-nucleus scattering 
based on the new formulation of the model should be the subject of a separate study, we have 
considered selectively in the present work the data of the elastic neutron scattering on the 
lead target, measured at the incident energies of 30.4 and 40 MeV 30]. Given no Coulomb 
interference and the energies close to the Fermi energy (Ep & 38.7 MeV obtained with the 
Fermi momentum kp 1.36 fm -1 ), the considered elastic neutron scattering data should 
be a good test ground for the improved single-folding approach suggested in the present 
work, with the rearrangement effects consistently taken into account. For the 208 Pb target 
we have used the empirical neutron and proton densities deduced from the high-precision 


elastic proton scattering at 800 MeV by Ray et al. 


60 


61j . These densities are available 
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FIG. 19. Local momentum k(E, R ) of the incident neutron (upper panel) and the scaling func¬ 
tion g(k(E, R)) (lower panel) determined self-consistently from the HF+RT real folded n+ 208 Pb 
potential obtained with the CDM3Y6 interaction at the neutron energy E = 30.4 MeV. 


the analytical form, and they have been used recently in our folding model analysis 62] of 
the charge-exchange ( 3 He,t) scattering to the isobar analog state of the target, to determine 
the thickness of the neutron skin in 208 Pb. 

From the folded n+ 208 Pb potentials shown in Fig. [17] one can see that the standard 
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folding method [19] gives a rather deep HF-type folded potential. After the RT is taken 
properly into account, the HF+RT folded potential becomes substantially shallower, much 


closer to the empirica 
of the nucleon OP 


Woods-Saxon potential given by the global CH89 parametrizations 
, which was proven to be accurate for the elastic nucleon scattering 


from medium and heavy targets at the incident energies below 100 MeV. The agreement of 
the HF+RT folded potential with the empirical CH89 potential becomes better, especially 
at the neutron energy of 40 MeV, when the scaling g(k(E } R)) function is consistently taken 
into account. To further test the folded n+ 208 Pb potentials shown in Fig. [TT] they were 
used (without any further renormalization of their strength) as the real OP in the standard 


optical model (OM) calculation of t 
the code ECIS06 written by Raynal 


re elastic n+ 208 Pb scattering at 30.4 and 40 MeV using 


taken from the global systematics CH89 
sections with the measured data 


63|, with the imaginary and spin-orbital parts of the OP 
From the comparison of the calculated cross 
one can see that a very good OM description 


30] in Fig. 


of the data has been obtained with the real folded OP after the RT is included. Although the 
OM fit is marginally improved after the HF+RT folded potential is scaled with the scaling 
function g(k(E, A)), it is important to note that this scaling function is resulted from the 
realistic momentum dependence of the s/p potential in the NM discussed in Sec. Q4] The 
local momentum k(E , R ) of the incident neutron is largest in the center, and approaches its 
asymptotic value at the potential surface (determined for a free neutron whose kinetic energy 
is equal the incident energy) (see upper panel of Fig. [T9]) . As a result, the g(k(E, R )) value 
is ranging smoothly from about 0.96 at small radii to unity at the potential surface, which 
is a rather mild mean-held effect on the shape of the folded nucleon-nucleus potential. In 
any case, such a scaling function is physically much more consistent compared to a constant 
g(E) factor used in the earlier version of the folding model [l9j at each nucleon incident 
energy. The standard HF folded potential is too deep in the center and cannot deliver a 
good OM description of the data (see FigJISli unless its strength is renormalized by a factor 
Nr fa 0.85. A renormalization factor Nr < 1 of the real folded OP was often obtained in 
the earlier folding model analyses of the elastic nucleon scattering 
contribution from the RT is likely the main reason. 


20 ], and the lack of the 


With a very good OM description of the considered elastic neutron scattering data by 
the HF+RT folded potential, the reliability of the folding model in predicting the nucleon- 
nucleus OP seems much improved. In view of the new folding formalism suggested in the 
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present work, a systematic folding model analysis of the elastic and inelastic nucleon-nucleus 
scattering over a wide range of energies should be of high interest. 


VI. SUMMARY 


A consistent HF study of the asymmetric NM has been done using the CDM3Y3 and 
CDM3Y6 density dependent versions of the M3Y-Paris interaction, with the focus on the 
rearrangement term of the s/p potential that arises naturally when the Hugenholtz-van Hove 
theorem is applied to calculate s/p energy from the total NM energy. Based on the exact 
expression of the RT of the isospin dependent s/p potential given by the HvH theorem at 
each NM density and the empirical energy dependence of the nucleon OP observed over a 
wide range of energies, a simple method has been proposed to account effectively for the 
momentum dependence of the RT of the s/p potential in the standard HF scheme, with an 
explicit contribution of the RT added to the density dependence of the CDM3Yn interaction 
(fTTTh The HF+RT s/p potentials obtained at the different NM densities and neutron-proton 


asymmetries agree reasonably with those predicted by the microscopic BH 
the NM that have the higher-order rearrangement term properly included 


calculations of 


411 - 


Given a direct link between the s/p potential in the NM and the nuclear symmetry 
4 7], we have determined anew the parameters of the isovector densiy dependence of the 
CDM3Yn interaction (fTTTl by matching our HF+RT results for the IV part of the nucleon 
OP in the NM with those of the BHF calculation by the JLM group 31j, [32]. Using these 
new parameters, the calculated nuclear symmetry energy S(p) agrees nicely with the latest 
empirical constraints at p < p 0 as well as the results of the ab-initio calculations of the 
asymmetric NM at p > /+• With the high-momentum part of the s/p potential modified 
based on the observed energy dependence of the nucleon OP, the nuclear symmetry energy 
S(p) obtained in the parabolic approximation from the difference between the neutron and 
proton s/p energies turns out to be in a very good agreement with that given by the ab- 
initio calculations over a wide range of the NM densities. This result indicates that one 
might indirectly learn about the density dependence of the nuclear symmetry energy from 
the extensive OM studies of the elastic nucleon-nucleus scattering at low, medium, and 
intermediate energies. 

The momentum dependence of the s/p potential obtained in the HF+RT calculation 
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has a very well defined linear 5-dependence of the n/p effective mass splitting (l33lh with 
m n-p(Poi $) ~ (0.26 ± 0.01)5 that agrees well with the empirical constraint from the recent 
analysis of the terrestrial nuclear physics experiments and astrophysical observations [dJ. 
The m*_ p (p, 5) value was found to be readily increased with the increasing NM density up 
to 2p 0 , while retaining its linear 5-dependence. Although there is no empirical constraint for 
the n/p effective mass splitting at the high NM densities to compare with, we have shown 
here that a proper treatment of the RT and a realistic momentum dependence of the single¬ 
nucleon potential are prerequisites for the determination of m* n _ p (p, 5) at the different NM 
densities that can be of interest for the nuclear astrophysical studies. 

A very important milestone of the present work is that the proper treatment of the 
rearrangement effects and momentum dependence of the nucleon mean-field potential in 
the HF calculation of the NM has lead us to the important physics inputs that enable a 
consistent inclusion of the RT into the HF-type folding model calculation of the nucleon 
OP of the finite nuclei in the same mean-field manner. The contribution of the RT has 
been shown, in an application of the extended folding model to study the elastic n+ 208 Pb 
scattering, to be vital in obtaining the realistic shape and strength of the real nucleon OP. 
The predicting power of the folding model for the nucleon OP seems much improved. A 
systematic folding model analysis of the elastic and inelastic nucleon-nucleus scattering over 
a wide range of energies is now planned with the extended folding formalism. 
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